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We demonstrate how to construct dynamic phase diagrams for complex fluids that undergo 
transitions under flow, in which the conserved composition variable and the broken-symmetry order 
parameter (nematic, smectic, crystalline, etc. ) are coupled to shear rate. Our construction relies on 
a selection criterion, the existence of a steady interface connecting two stable homogeneous states. 
We use the (generalized) Doi model of lyotropic nematic liquid crystals as a model system, but the 
method can be easily applied to other systems, provided non-local effects are included. 



Complex fluids in shear flow display a range of be- 
haviors which is only beginning to be unearthed [p}-p|. 
Shear can perturb equihbrium phase transitions (e.g. 
the isotropic-to-nematic (I-N) liquid crystalhne 
and isotropic-to-lamellar transitions), and induce 
structures, (e.g. bilayer onions that exist only as 
metastable equilibrium phases. A related phenomenon is 
dynamic instability in non-Newtonian fluids, when the 
theoretical homogeneous stress-strain-rate constitutive 
relation exhibits multi-valued behavior, as in theories 
of polymer melts |ll|,|l^ and worm- hke micelles 
Such models are often used to describe, for example, the 
spurt efl'ect, whereby the flow rate of a fluid in a pipe 
changes discontinuously as a function of applied pressure 
drop . The most important unresolved question about 
non-monotonic flow curves, such as those in Fig. |l|a, is: 
what determines the stress, if any, at which the system 
phase separates into 'bands' ? Suggestions have included 
(a) variational hypotheses (Ref. (b) assuming the 

stress at the top of the stable viscous branch is selected 
("top jump") (c) including geometrical effects 

| fi8[ ; and (d) incorporating (physically present) non-local 
contributions to the stress [^| Jl9|j20[| . In this Letter we 
pursue (d) and explore in detail the utility of constitu- 
tive inhomogeneous effects in resolving the issue of stress 
selection in complex fluids, which has occupied the rhe- 
ology and physics communities in the guise of either un- 
stable flows [0 or non-equilibrium phase transitions [D . 

After a general discussion we introduce semi-phenom- 
enological equations of motion for rod-like molecules in 
solution, extending the Doi model to inhomogeneous 
flows. Our study has the following goals and results: (1) 
We present a general recipe for computing phase separa- 
tion under flow, and hence the experimentally measured 
rhcological behavior; (2) We point out that the proper 
field variable {either stress or strain rate) may not be 
unique, a feature absent from equilibrium systems; (3) 
Using concepts developed for dynamical systems theory 
we conclude that stress selection of models with non- 
local (in space) differential constitutive relations in pla- 
nar shear flow is unique; i.e., as with equilibrium phase 
transitions, it occurs along a hyper-surface of lower di- 



mension than that of the field variable space. Our dis- 
cussion is facilitated by examining the stress-strain-rate- 
composition surface, a representation we have not seen 
before and hope will become commonplace. 
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FIG. 1. (a) Ifomogeneous stress ax 



strain rate 7 be- 



havior for various excluded volumes (u = cj)Lv2/ct) as cal- 
culated from the Doi model with quadratic closure [21,23]. 
We use the dimensionless parameters 7 = / {dDroViv^) a-nd 
Bxy = (7xyV2iog L/{9nrisDro)- Dotted lines mark unstable 
branches. The plane in (b) is at fixed stress. At Ue both 
isotropic and nematic phases are metastable in zero shear. 
Here and below L — 5. 



The curves in Fig. ^ are reminiscent of pressure- 
volume isotherms in a liquid-gas system and, accord- 
ingly, we seek to construct a 'phase diagram' by pur- 
suing an analogy between homogeneous stable steady 
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states and equilibrium phases. As in equilibrium, non- 
equilibrium 'phases' may be separated, in field variable 
space, by hyper-surfaces representing continuous (e.g. 
critical points) or discontinuous ('first-order') transitions 
1^,0. Coexistence implies an inhomogeneous state span- 
ning separate branches of the homogeneous fiow curves. 
In contrast to equilibrium bulk phases, where field vari- 
ables are uniquely identified, here we may consider phase 
coexistence at a given shear stress, with the interface 
parallel to the vorticity- velocity plane (curves b—e); or a 
given strain rate (curves c— e), with the interface parallel 
to the velocity-velocity gradient plane. The appropri- 
ate field variable is thus determined by the nature of the 
constitutive relation (for curve b, common strain rate is 
impossible) , or even by the flow history or the rheometer 
(e.g. fixing the stress or strain rate in curves c— e) . In this 
Letter we compute coexistence at a common stress, but 
note that computation for a common strain rate is anal- 
ogous. In equilibrium, minimizing the total free energy 
leads to equality of field variables between two phases and 
the Maxwell common tangent condition (e.g. the equal 
osmotic pressure condition, aided by equal chemical po- 
tential, in rod suspensions p^). In dynamics, the former 
follows from the equations of motion plus the stationary 
condition, with an extra shear degree of freedom. How- 
ever, the lack of a criterion to replace the latter leaves 
an unresolved degeneracy (e.g. which stress is selected 
for a given averaged strain rate). Flow experiments on 
worm-like micelles find a well-defined transition stress for 
the onset of banded flows |||,§,||,D, which suggests that 
the degeneracy is not physical. 

For the thermotropic (i.e. restricted to the melt com- 
position) I-N transition under flow, Olmsted and Gold- 
bart resolved this degeneracy by rejecting those pos- 
sibilities which did not admit a stationary interface solu- 
tion and found, numerically, an apparently uniquely se- 
lected state. They included gradient terms, which penal- 
ize (energetically) spatial microstructure variations and 
dictate the interface structure. The importance of gra- 
dient terms was also recognized in Refs. p3| , p^ 20 1. The 
existence condition of a stationary interface selects one 
among a band of possible coexisting solutions, and is fun- 
damentally different from augmenting a local constitutive 



161 or assuming selec- 
17|. Below we justify 



model with a variational principle 
tion at the limits of stability | p2| , p^ f 
the uniqueness of a selected state. 

Now we proceed with our calculation. The free energy 
of a solution of rod-like molecules of volume fraction (j) in 
an athermal solvent (e.g. as in Ref. pS]) is given by 



.F(0, Q)- fc«T / <j ^ log0 log (1^ 



i (l-iu)TrQ2_ luTiQ 
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+ ^K{S/^Qp^f]+y{Vcj,fY (1) 
Here Vr and are rod and solvent monomer volumes; 



Qap is the nematic order parameter tensor ||2^; and 
u = 4'Lv2/a is Doi's excluded volume parameter [^ , 
where L is the rod aspect ratio and V2 and a are 0(1) 
geometrical prefactors. This free energy includes the en- 
tropy of mixing, the oricntational contribution yielding 
an I-N transition |^l[ , and spatial correlations of compo- 
sition and nematic order within the one constant approx- 
imation. The phenomenological gradient terms may, in 
principle, be calculated from a microscopic model. 



The equations of motion are 



microscopi 

imii: 
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with Kajj — V/3Wq and p the density. The stress tensor is 

cr + 2?/K^ + cr^et,(0,Q), (5) 

where = (k^^ -|- K/3a)/2 and we take 77 to be the sol- 
vent viscosity rfg for simplicity. The reversible stress due 
to the nematic order is pl|,p|,p4[ 



m + H Q Q H VQc^p 



(6) 



where H ~ —5J-/5Q. The isotropic stress from the force 
density {^4))5T / Scf) |Q has been neglected. In Eq. (||) 
the (reactive) ordering term F is given by 

F{k.,Q)^Ik' + k-Q+Q-h'^ -2{Q+\l)Tr{Q-K). (7) 

For simplicity, we have chosen the form appropriate for 
an infinite aspect ratio |21|. The dissipative portion G is 



«(</>, Q) 



(8) 



where D^o is the single-rod rotational diffusion coeffi- 
cient, vi is a geometric prefactor, and t is the rod length 
|^l|. The chemical potential /i drives the current J, 



J = M ■ V/i, 



(9) 



where M is the mobility tensor and p = 5J- / 5(j}. 

For other systems, equations like Eqs. (^) and (^ gov- 
ern the conserved and broken-symmetry (or other long- 
lived) variables, respectively. For some local models, in- 
ternal dynamics (Eq. ^ can be eliminated to give the 
stress as a history integral over the strain rate. In poly- 
mer melts jll|] , and worm- like micelles 1 14 far from a ne- 
matic regime, this leads to non-monotonic stress-strain- 
rate curves. 

We seek stable steady-state solutions to Eqs. (||-^) for 
planar shear, v(r) = 72/x. Integrating Eqs. (||,^) along y 
yields axy{4>,j,Q) — <Jo and p,{(j),Q) — /io, where do is 
the applied stress. One integration constant of Eq. (^) 
is zero from the boundary condition Jy = 0, while /io is 
determined below P7|. For homogeneous solutions, Q 
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may be eliminated from a^y and /i using Eq. (|). The 
stress is shown in Fig. |^. Because !F{(p, Q) describes an 
I-N transition, multiple roots for Q may exist at a given 
stress MM, with distinct strain rates. 
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FIG. 2. (a) Reduced strain rate 7(11) for the stress contour 
in Fig. (b) fj-{u) along the curve in (a). A selected pair 
has been indicated in (a) and (b). The inset shows another 
(topologically) possible fi{u). Dashed lines mark the limits of 
stability of coexisting pairs. 



ports the existence of such a solution in the modified Doi 
model, and is stable against perturbations in y. Because 
non-local effects (see Ref. ||2^,|l^ for diffusion effects) al- 
ways exist in reality, and pathological degeneracies only 
occur in local models, we expect that models that can 
resolve the interface structure have unambiguous phase 
diagrams. 
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FIG. 3. Phase diagram in the (axy—u) (a) and (7— u) (b 
planes. Broken lines trace the loci of points 1 — 4 in Fig, 
and the solid lines are calculated tie lines. 



For coexistence at fixed stress ctq (the interface has a 
normal parallel to y), a stress contour traces a line in the 
(7 — 0) plane (Fig 0a). We plot /i(0) along this line in 
Fig. yb. Coexisting phases must have the same /i. For 
concentrations (j) outside points 1 and 4, /iq is determined 
uniquely. Between these points, we must determine /iq 
at which two phases coexist. Following Ref. we select 
the fiQ which permits a stationary inhomogeneous solu- 
tion satisfying (7j;y = ao, ^ = ^0 f^nd Eq.(|3|) {0 = F + G). 
This is not an auxiliary assumption, but follows from the 
inhomogeneous equations of motion. Given a selected 
/io and a mean concentration (j), the portions of the two 
phases are fixed by the lever rule, and the mean shear 
rate can be calculated. In the static limit (7 — 0) the 
stationary interface condition follows from minimizing 
!F{4>, Q) with a given (j) constraint. Functional minimiza- 
tion with respect to (/)(?/) allows the interface position to 
move, thereby adjusting the ratio of the two phases and 
recovering the common tangent construction. 

To see how ^0 can be selected uniquely, note that the 
stationary solutions to Eq.(||), axy = ctq, and fi = fiQ, 
are nonlinear ordinary differential equations (ODE) in- 
volving d/dy. These ODE may be converted into an 
equivalent set of first order ODE with dependent vari- 
ables {(t>,d(j)/dy,j,Q,dQ/dy}. In the ODE phase space, 
the interfacial solution corresponds to a trajectory (a 
'heteroclinic saddle connection') joining two fixed points 
(the homogeneous states). As /io changes for fixed co, 
the phase flow changes catastrophically at those isolated 
values of fio where the desired trajectory exists. [This 
is proven for differential constitutive relations in planar 
flow by showing that the saddle connection, assum- 
ing it exists, is of the non-transverse type [^]. This ex- 
plains the apparent uniqueness of a selected /io, given ao, 
found numerically. Our (one-dimensional) solution sup- 




The selection criterion determines the tie lines on 
and varying the stress yields the phase diagram of 
|33 ||2^. In the latter, for 7 = the tie line is horizon- 
tal, and for 7 > tie lines have positive slope, because 
the more concentrated nematic phase has a lower effec- 
tive viscosity. For models where fi{(f>) has the shape of 
the inset of Fig. ||b, tie lines have negative slopes. 
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FIG. 4. axy vs. 7 for various compositions. To read this 
from Fig. ^, choose a u and increase the stress. In the 
two-phase region jump from stress tie line to stress tie line 
until the one phase region is reached. 

Fig. H shows the stress-averaged-strain-rate curves as 
would be measured in an experiment. Here the co- 
existence region comprises lines with positive slopes. 
For models with vertical tie lines in Fig. indicat- 
ing a composition-independent transition (as in the Doi- 
Edwards model of polymer melts), the plateau would be 
flat. This is a graphical expression of the explanation of 
a sloped plateau given by Schmitt et al. [|7| (however, 
they assumed "top-jumping"). In general one must vary 
both average concentration and strain rate to maintain 
constant shear stress, a situation familiar from equilib- 
rium multi-component phase coexistence. The 'plateau' 
need not be a straight line; its shape is determined by the 
change in splay of the tie lines with increasing strain rate 
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in Fig. Experiments on wormlike micelles g] have 
found the plateaus to be steady states and the 'spines' to 
be metastable branches. 

In summary, the steps to compute phase separation un- 
der flow are: (1) Determine the full inhomogeneous equa- 
tions of motion; (2) Analyze stable homogeneous flows. 

(3) Choose the candidate field variable for phase coexis- 
tence. Multiple choices (e.g. stress or strain rate, accord- 
ing to interface geometry) must be analyzed separately. 

(4) Identify pairs of phases with the same field variables, 
as in Fig. ^ and Fig. |^b for axy and fj,. (5) Determine the 
phase diagram (in field variables) by requiring a station- 
ary interface between homogeneous solutions. To com- 
pare with experiments, density variables can be used, as 
in Fig. ^, to find (6) the tie lines where the lever rule 
applies and (7) the space-averaged stress-strain-rate re- 
lations (Fig. U). Other multi-component systems retain 
the structure of Eqs. (|[]^), with Eq. (^ governing the 
dynamics of some structural order parameter. 

We have calculated coexistence at fixed stress, while 
curves c—e in Fig. |l]a admit, in principle, coexistence at 
fixed strain rate. A full solution of Eqs. requires 
analyzing both interface directions, which is reminiscent 
of, and may be relevant to, the problem of the orienta- 
tion of diblock lamellar phases in shear |^^. Also, we 
have only considered flow-aligning solutions to the Doi 
model. Another, so-called 'log-rolling', solution is sta- 
ble at high shear rates (29) and will be included in a 
complete treatment [Q. The extended Doi model is ap- 
propriate for hard rod suspensions |^^ , and we welcome 
experiments on these systems. The flow instability is due 
to perturbation of an equilibrium phase transition; while 
systems such as worm-like micelles probably possess some 
combination of perturbed (I-N) and dynamic transitions 
P,p^ , which could yield a stress-strain-rate-composition 
surface (Fig. |l|b) with multiple folds. 

We have not considered the important problem of the 
stability of an undulating interface, which may restrict 
the choice of field variables ||3^]. Stability analysis by 
Renardy jsT) on banded flow in the Johnson-Segelmann 
model (which has the qualitative behavior of curve b in 
Fig. |l|a, but no gradient terms) revealed a stationary 
interface at any stress in the two-stress region, and a 
family of unstable high-wavenumber undulations. Gra- 
dient terms break the stress degeneracy and 
introduce a stabilizing surface tension and dissipation 
within the interface as it moves, which affects the stabil- 
ity analysis. Finally, complementary to the planar shear 
case analyzed above, Greco and Ball ||l^, have recently 
demonstrated the important result that, for a Johnson- 
Segelmann fluid in cylindrical Couette flow, a stationary 
interface exists at a selected stress: coexistence is appar- 
ently influenced by the curved boundary geometry of the 
flow. 
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